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A simplified model of particle transport at a quasiparal- 
lel one-dimensional collisionless shock is suggested. In this 
model the MHD-turbulence behind the shock is dominated 
by a circularly polarized, large amplitude Alfven wave origi- 
nated upstream from the turbulence excited by particles leak- 
ing from the downstream medium. It is argued that such a 
wave having significantly increased its magnetic field during 
the transmission through the shock interface can effectively 
trap thermal ions, regulating their leakage upstream. To- 
gether with a background turbulence this wave also plays a 
fundamental role in thermalization of the incoming ion flow. 
The spectrum of leaking particles and the amplitude of the 
wave excited by these particles are selfconsistently calculated. 
The injection rate into the first order Fermi acceleration based 
on this leakage mechanism is obtained and compared with 
computer simulations. The related problem of shock energy 
distribution between thermal and nonthermal components of 
the shocked plasma is discussed. The chemical composition 
of the leaking particles is studied. 



I. INTRODUCTION 

The problem of energy dissipation in collisionless 
shocks in plasmas is old and exceedingly difficult 
Moreover, there exist persuasive theoretical arguments 
p0| , pl| corroborated by numerous simulations ( e.g., |12|) 
that this problem cannot be solved by considering exclu- 
sively the thermalization of the bulk plasma flow when 
the latter passes through the shock. A significant part if 
not almost all the energy of a strong large shock may be 
channeled into a small minority of particles accelerated 
through multiple crossing of its front. This acceleration 
mechanism is known as the first order Fermi or diffusive 
shock acceleration. However, this does not circumvent 
the problem of collisionless thermalization. The reason 
is that a small fraction of thermal ions that leak or re- 
flect from the shock play a crucial role in the collisionless 
energy exchange between the bulk upstream motion and 
thermal and/or nonthermal (accelerated) components of 
the downstream plasma. These ions generate waves in 
the foreshock region whose growthrate and amplitude 
arc directly related to their density. Yet they provide a 
seed, or injection population for the further acceleration. 
These two aspects of the shock dissipation are clearly in- 
terrelated. By wave excitation these particles create a 
scattering environment allowing them to cross the shock 



repeatedly which is necessary for the Fermi mechanism 
to work. 

One of the most important parameters of collisionless 
shocks is the angle n B between an ambient magnetic 
field and the shock normal. While so called perpendicular 
shocks (0„b — tt/2) should clearly have a distinct shock 
transition because the hot downstream plasma cannot 
penetrate upstream for more than one ion gyroradius, 
their parallel counterparts (8 n B *C 1) are not so suitable 
for confinement of the heated downstream plasma since 
it may penetrate far upstream moving along the field 
lines. We will confine our consideration below to this 
latter category of collisionless shocks. In general, shocks 
with n B < 7r/4 are somewhat superficially referred to 
as quasiparallel whereas the rest (tt/4 < 0„b < tt/2) as 
quasi-perpendicular . 

It is important to realize that the distribution func- 
tion of the backstreaming particles cannot be inferred 
solely from the macroscopic parameters of the down- 
stream plasma even if the thermalization mechanism is 
properly understood. One obvious reason for this is the 
following. Among the backstreaming particles one can 
find not only those which simply arose from the random- 
ization of the upstream flow at the shock (as backscat- 
tered from the downstream medium or reflected from the 
shock interface) but also the particles of these two types 
which crossed the shock more than once and therefore 
have gained some energy fl3| . 

The above arguments suggest that the injection prob- 
lem can be divided into the following two tasks: (i) given 
the shock conditions one determines the distribution of 
particles originating from the upstream flow after they 
crossed the shock for the first time. Subsequently, one 
identifies those particles which are capable of crossing the 
shock also in reverse direction (first generation of injected 
particles), (ii) One follows the (stochastic) trajectories of 
these particles when they multiply recross the shock un- 
til they are swept downstream or have achieved energies 
acceptable for the standard description of diffusive shock 
acceleration (see e.g., Ref. jl4j or |L5j for a review). 
The first task belongs to collisionless shock physics 
and can at least formally be treated independently 
of the diffusive shock acceleration process. The second 
constitutes the injection problem itself as a part of dif- 
fusive shock acceleration theory and can be formulated 
in more detail as follows. Suppose the task (i) is solved. 
Then, given the distribution of thermal particles that are 
able to penetrate into the upstream region, one calculates 
the high energy asymptotics of their distribution. This 
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provides the coefficient in the power-law solution of the 
standard acceleration theory and thus the injection rate. 
It cannot be obtained within the standard theory since 
the latter is unable to describe low energy particles with 
anisotropic pitch angle distribution. 

The solution of the injection problem (ii) as formulated 
above has been obtained analytically in Ref. |13|. The 
high energy asymptotics of this solution indeed matches 
the power-law of the standard theory. At the lower en- 
ergy end it smoothly joins the downstream thermal dis- 
tribution. This thermal distribution that determines the 
flux of leaking particles has been simply assumed to be 
created by a unspecified shock randomization process, in 
other words due to the action of a "thermostat" . In this 
paper we suggest a simplified thermostat model. In the 
next subsection we put it in the general context of colli- 
sionless shocks. 



A. Thermostat model 

Some of the thermostat properties are readily known 
from ordinary Rankine-Hugoniot (RH) jump conditions 
that constitute the conservation of mass, momentum and 
energy fluxes across the shock jl6| . In the simplest case 
when both the magnetic field and the flow velocity are 
perpendicular to the shock front the RH conditions are 
the same as in ordinary gasdynamics 



beam would be trapped in this wave and carried back 
downstream. Therefore, considering strong shocks with 



+ Pl 



P2U2 
P2U2 



P2 



1 3 

2^ + 7" 



1 



7-1 



P2U2 



(1) 
(2) 

(3) 



Here the index 1 (2) refers to the upstream (downstream) 
medium, p ~ m p n, u, and P are the mass density, ve- 
locity and the gas kinetic pressure, respectively, m p is 
the proton mass. The adiabatic index 7 can be set to 
7 = 5/3. The pressure and energy of the MHD waves 
that are implied to participate in the shock thermal- 
ization process as a substitute for the binary collisions 
are neglected here for the following reasons. First, if 
the initial state (1) refers to the far upstream region in 
which the backscattering particles do not penetrate and 
the waves are thus not present there, the left hand sides 
of the above equations are written correctly. If the final 
state (2) is also taken so far downstream where the ther- 
malization process is completed and the energy of the 
waves excited at the shock interface is damped, the same 
is true for the right hand sides. Even if the final and 
initial states are taken substantially closer to the shock 
interface these equations correctly describe the jump re- 
lations in the case of a very strong shock. Indeed, the 
waves that are excited in a highly supersonic and su- 
peralfvenic upstream flow via cyclotron resonance with 
the beam of backstreaming particles, are the transverse 
(SB _L Bo) MHD waves and their amplitude SB can- 
not exceed the unperturbed field B since otherwise the 



u\/v\ = A-Kp\u\/Bl > 1 we may neglect the 



contribution of the wave energy upstream. Downstream, 
it is increased by a factor ~ {U1/U2) 2 due to the shock 
compression of the perpendicular component of the wave 
magnetic field but the wave energy may still be neglected 
compared to the thermal and the bulk motion energy in 
strong shocks with Ma 3> 1. For whatever reasons the 
RH conditions have been observationally proven to be 
quite reliable in collisionless shocks |lj||l8| ■ 

In the case of high Mach number shocks A/J = 
u\p\/"fP\ ^> 1 the RH conditions simply yield (for 
7 = 5/3): u\ju2 = 4 and the downstream temperature 
T2 = im v u\. If one is assuming that all the downstream 
particles that have velocities v z < —U2 (where v z denotes 
the normal to the shock front velocity component, and 
the shock is moving in negative z— direction) can cross 
the shock, this information about the thermostat suffices 
to obtain an injection solution. On the other hand, such 
an assumption can only give an upper bound on the in- 
jection efficiency and it highly overestimates the latter. 
In principle, many particles in such a high temperature 
downstream distribution have negative velocities in the 
shock frame of reference and, therefore, can potentially 
escape upstream. But this would be a purely kinematic 
picture. In reality too intensive escape would result in a 
fast excitation of waves scattering the beam back down- 
stream over a short distance. This process can be un- 
derstood within the framework of the quasilinear theory 
of the cyclotron instability developed for hom ogeneous 
plasma in Ref. JTi| (see also ]l3| and Sec. VII below). 

Consider the phase space of such a beam emanating 
from the shock upstream on the plane (v\\,v±) in the 
shock frame of reference Fig. [l], where , are paral- 
lel and perpendicular to the shock normal (and to the un- 
perturbed magnetic field) components of particle veloc- 
ity. The arrangement of the beam in the phase space cor- 
responds qualitatively to its parameters calculated later 
in this paper. As soon as the beam appears on the up- 
stream side of the shock it starts to generate MHD waves 
that move to the left in the local plasma frame at the 
Alfvcn speed — Ca but are in fact continuously convected 
back to the shock with the main flow since m ^> Ca- 
Their resonance length is of the order of ui/w c i where 
uj c i is the ion cyclotron frequency. The beam itself is be- 
ing scattered in pitch angle by this selfgenerated waves 
around the center at u\ — Ca so that it gradually dif- 
fuse into the Um > region which means that it returns 
downstream. The relaxation length of the beam, i.e., the 
depth of its penetration upstream may be estimated as 
Ir ~ u i/l where 7 is the growth rate of the cyclotron in- 
stability upstream. Estimating 7 from quasilinear equa- 
tions ( e.g., Ref. one may obtain for Zr, 
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where lu p i is the ion plasma frequency, n^, hq are the 
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beam and the background densities, respectively, and 
Avu is the beam width in vu . The last factor in Zr is 
in fact only slightly less than unity so that the upstream 
part of the entire shock transition is roughly no/^b > 1 
ion inertial lengths, c/uj v \. 

In principle the turbulence generated by the beam and 
growing in the downstream direction could gradually sat- 
urate at a distance > Ir and the plasma flow beyond this 
distance could be declared as a downstream state. How- 
ever, this is not what actually happens as both the obser- 
vations and numerical simulations of strong shocks reveal 
(see e.g., 0, ||). Namely, there exists a relatively sharp 
shock transition inside of this structure where the am- 
plitude of magnetic pulsations increases over a distance 
~ c/a->pi, the bulk flow slows down nearly to its down- 
stream value and particle orbits spread in velocity to an 
approximately thermal width in a substantially increased 
magnetic field. 

Existing ID and 2D hybrid simulations (ions treated as 
particles and electrons as a fluid) of quasiparallel shocks 
H^O - 24 1 also provide some clue to how particles are con- 
fined on the downstream side of the shock. Namely, 
despite the fact that there are indeed many particles 
just downstream of the shock front that have velocities 
«ll < 0, only few of them recross the shock. This suggests 
that the heated downstream plasma is directly locked by 
the waves which have been excited in the upstream region 
and transmitted then across the shock. In part, they may 
be generated at the shock interface. Because of a sub- 
stantial downstream increase of the wave amplitude due 
to the shock compression, the perpendicular component 
of B— field may become large enough to trap the bulk of 
the particles and to sweep them downstream. An addi- 
tional factor that can also help to confine the shocked 
plasma is an electrostatic barrier that is also observed 
in simulations. The analysis of the backscattered ions 
performed by Quest || shows that the upstream flux of 
these particles is significantly smaller than it would be if 
all downstream particles with negative velocities were to 
stream freely across the shock. 

Another important conclusion that can be drawn from 
the ID hybrid simulations (see e.g., Ref. is that the 
downstream wave turbulence is dominated by a mode 
excited by the backstreaming beam in the upstream re- 
gion. One can think of a circularly polarized quasi- 
monochromatic MHD wave convected downstream with 
the bulk plasma and extended there over a distance of 
the order of hundred ion inertial lengths. The amplitude 
of this wave is quite large SB /Bo ^3 — 4 (see also Ref. 
p3| ) . Traces of the quasi-regular particle motion in the 
wave are also seen in the simulations Q. Beyond this 
distance the wave is gradually damped and the particle 
distribution tends to a thermal one. However, unlike the 
upstream part of the total shock transition the dis- 
tance where it happens is very difficult to assess analyt- 
ically. The reason is that the level of turbulence is very 
high and a weakly turbulent approach is probably im- 
possible (see, however, Sec. III). This part of the shock 



transition should rather be regarded as a BGK type of 
wave |p5| . Such an attempt is made in a companion paper 
p6[ where we consider the slowing down and the heating 
of the upstream plasma caused by interaction with this 
wave downstream. 

The present paper suggests a model of ion leakage al- 
lowing the determination of the thermostat production 
rate that was only parameterized in the injection the- 
ory [Fq ]. The starting point of the model is the parti- 
cle dynamics in the Alfven wave behind the shock p7| . 
The particle phase space is divided into two parts. One 
part contains the "trapped" particles which are convected 
downstream with the wave. The second part contains 
"untrapped" or adiabatically with the wave interacting 
particles, and particles trapped in the nonlinear Doppler 
resonance. These particles can escape upstream when 
their averaged velocity with respect to the wave is high 
enough. 

It is not assumed, however, that particle motion in 
this model is purely regular. Instead, a downstream tur- 
bulence that should exist together with the monochro- 
matic wave allows particles to cross the boundary be- 
tween trapped and untrapped regions and also between 
the "lock" and "escape" states of particles in the phase 
space. Although this motion across the invariant mani- 
folds of the regular dynamics is assumed to be relatively 
slow it should produce the necessary entropy by virtue of 
the standard arguments i t^lpsfl : even a very weak pertur- 
bation can quickly randomize the particle motion while 
imposed onto the regular motion in the large amplitude 
wave since the main job is done by this wave via fast 
phase mixing. From this reasoning we obtain the fraction 
of the downstream particles that escape upstream from 
such a trap as a function of their energy. This allows 
us to calculate the unknown "thermostat" distribution 
function in terms of the wave amplitude. 

The final step of this scheme is the determination of the 
wave amplitude from the beam density. This provides, 
in fact, an equation for the wave amplitude since the 
beam density depends on this as well. Upon solution of 
this equation one obtains both the beam density and the 
wave amplitude as functions of shock parameters. 

In the next section we briefly review the particle dy- 
namics in a monochromatic circularly polarized wave, in- 
troduce suitable variables, and consider the special case 
of a very strong wave (SB/Bq ^> 1). In Sec. [n| the role of 
the background turbulence is discussed and adiabatically 
leaking particles are identified. In Sec. fv\ the probability 
of their escape is calculated as a function of energy and 
wave amplitude. In Sec. |v| the resonant escape is consid- 
ered. Sec. VI deals with the dependenc e of escape fluxes 
upon the mass to charge ratio. In Sec. VII we calculate 
the wave amplitude upstream. This enables us to obtain 
the particle flux without parameterization, simply as a 
function of the Mach number. After a bri ef su mmary of 
the considered leakage mechanism in Sec. VIII we calcu- 
late in Sec. IX the spectrum injected into the first order 



Fermi acceleration and compare it with the hybrid sim- 
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ulations. We conclude with a discussion of alternative 
mechanisms and possible applications to the problem of 
energy partition of collisionless shocks between thermal 
and nonthermal particles. 



II. PARTICLE DYNAMICS IN THE WAVE 

The equations of motion in the case considered here 
are fully equivalent to the equations of particle motion 
in a whistler wave. These equations have been exten- 
sively studied in the past. It is well known ^] that 
they are completely integrable. However, unlike the situ- 
ation with the whistler wave where one usually assumes 
SB /B < 1, we must concentrate on the opposite case 
as it was described in the preceding section. We there- 
fore present here a technically different description of the 
particle dynamics which is suitable to our purposes. 

Let us assume that the wave propagates in z- direction, 
i.e. normal to the shock and Bq = B z = const. We 
represent the total magnetic field in the reference frame 
moving with the wave in the following form 

B = B, z e z + B± (— e x cos koZ + e y sin k^z) (4) 

where {e;} denotes the standard basis in coordinate space 
and fco is the wave number. The electric field vanishes in 
the wave frame of reference and the equations of motion 
read 



dv e _ 

— = v x B 

at rripC 

dz 

— = v 7 
dt 



(5) 
(6) 



where the particle velocity v = v x e x + v y e y + v z e z . In- 
terested in the case 



B z 1 
£ ^Bl <<1 ' 

we introduce the cyclotron frequency 

eB ± 

uj ± = 

rripC 

and rescale the variables in Eqs. (pUq) as follows 



k z — > z, uj±t — > t, 



fc v 



Using these new variables Eq.(|5|) rewrites as: 
dv 



dt 



= v x (— e x cos z + e y sin z + ee z 



(J) 



(8) 



(9) 



(10) 



and Eq. (J6J) remains unchanged. It is convenient 
to make a further transformation in these equations 
(v x ,v y ,v z ,z) — > (X,q,v z ,z), in which 



It is easy to see that A is conserved, dX/dt = 0, and we 
arrive at the following one dimensional (i.e. integrable) 
dynamical system in the variables (q, v z ) 



dq e cos q 

dt~ Vz+ VA 2 - 2ev z ~ £ 

dVy, 



dt 



-\/X 2 — 2ev z si 



sing 



(12) 



The integral A 2 can be written in the old variables as 



A = 



2v x cos z + 2v Y sin z + 2ev z + 1 



Besides A there is another obvious integral, the energy 
v 2 — v 2 +v 2 +v z . Note that A 2 in Eq.(|l2|) can be negative 
in certain parts of the phase space; positive definite is 
only the quantity A 2 — 2ev z . However, being interested in 
the case e -C 1, we start our consideration of the system 
( |l2l ) from the simplest situation where A 2 3> \sv z \. As 
we will see the escaping particles interact with the wave 
adiabatically (v > 1) in this case. We shall return to the 
case of small and negative A 2 < e (resonant escape) in 
Sec. 0. 



A. Adiabatic wave-particle interaction 



Introducing the new variable 



V = v% + - cos q — £ 
A 



(13) 



and retaining only the terms of zeroth and first order in 
e, the system (Q) describes a simple pendulum 



dq drj 

-dt =m dt = - Xsmq 



with the Hamiltonian 



ff = 2Asin 2 (f) + ^ 2 . 



(14) 



(15) 



that is connected with the integrals v and A by means of 
the following relation v 2 = 2H+(X— l) 2 . It is convenient 
to introduce the standard action-angle variables in the 
system (II30), using the truncated action 



S = I rjdq 



v x + iv y = — \J X 2 — 2ev z exp [i(q — z)\ + cxp(-iz) (11) 



(16) 



as a generating function. The function k 2 = 2X/H = 
4A/[w 2 — (A — l) 2 ] divides particle phase space into two 
parts that will be superficially referred to as the region 
of trapped (k > 1) and untrapped (k < 1) particles. 
It should be pointed out that the particles with k > 1 
are not really trapped in the usual sense because they 
can become untrapped (k < 1) without changing their 
energy v 2 . We shall return to this point below. Here we 
only note that the analogy to e.g., the particle dynamics 
in a monochromatic Langmuir wave is incomplete in 
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this respect. In many other respects "trapped" particles 
behave as such, in particular their averaged velocity fj is 
zero. According to Eq. ( [i"3| ) this means that v z ~ e. 

For the untrapped particles, using equations ( |l5| , ^6|) , 
we obtain 



5 = 4 



x*(§.*) 



(17) 



where E is the incomplete elliptic integral of the second 
kind. It is convenient to define an action J as 



— sign(?7) 
27T JO 



rjdq 



so that we finally obtain 



J = -^sign(77)E(fc), 
7rfc 



(18) 



(19) 



where E denotes the complete elliptic integral of the sec- 
ond kind. Thus, the untrapped particles occupy the re- 
gions |J| > Sy/X/ir = Js, and far from the separatrix 
(k = 1) where k — > one simply has J — * rj. The angle 
variable tp conjugate to J is 



4, 



dS_ 



K(fc) 



(20) 



where F and K are the incomplete and complete elliptic 
integrals of the first kind. Note that tp — » q as k — > 0. 

We will not use the corresponding action-angle vari- 
ables for trapped particles. 



III. THE BACKGROUND TURBULENCE AND 
ESCAPE UPSTREAM 

The simple particle dynamics in a monochromatic 
wave considered in the preceding section implies that the 
amplitude of this wave is constant in space and time. 
This is certainly not the case for the wave associated 
with a shock. As we mentioned in the Introduction the 
wave is extended over a finite distance on the downstream 
side of the shock (we assume between z = and z = L, 
and k^L 3* 1) and decays at larger z. Thus, particle 
interaction with this wave occurs in the following way. 
First, particles that cross z = plane from z < become 
trapped (at least an appreciable part of them) and move 
downstream. Indeed, after crossing z — they feel a 
strong quasi-perpendicular wave field. The wavenumber 
can be estimated as 



U2 



ku 



(21) 



where fc u is the wave number of the most unstable and 
presumably strongest mode in the upstream region ex- 
cited by the escaping beam due to the cyclotron res- 



variables). As we shall see, uy in this resonant condi- 
tion may noticeably exceed u\. Here, we estimate k u as 
k u ^ wm/ui- Thus, the gyroradius of particles crossing 
the shock downstream is smaller than the wave length 



ko 



m - u 2 



< 



Ui - u 2 

U 2 



e < 1 



(22) 



onance kv 



= (we use here unnormalized 



and these particles must be effectively deflected in the 
wave magnetic field. Since the wave is not really 
monochromatic and low amplitude turbulence is always 
present as well, the motion of particles is not fully deter- 
ministic. It should be noted here that the background 
turbulence can be very well generated by particles them- 
selves. One generation scenario (the so called side band 
instability) has been studied for Langmuir waves in Ref. 
J32|, and for whistler waves in Ref. |33|. The results 
can be summarized briefly as follows. After the wave is 
switched on (or, what is more appropriate for our case, 
the plasma has entered the region of wave localization) 
and the particles have bounced a few times in the wave 
field, their distribution becomes "ergodic" and depends 
only on the action J. This happens due to the fast mix- 
ing in the phase variable if). This ergodic distribution 
is, however, usually unstable with respect to the excita- 
tion of satellites of the main wave due to the resonances 
u>k = nQ(J). Here, Uk and k ~ ko are the frequency and 
the wave number of a satellite in the main wave frame, 
f2(J) = dH/dJ is the frequency of particle oscillations in 
the wave and n is an integral number. 

A quasilinear theory of the backreaction of excited 
satellite turbulence on the main wave and on the ergodic 
particle distribution has been developed by this author 
p4j . This theory shows that the satellites change sig- 
nificantly both the particle distribution and the wave. 
However, they evolve relatively slowly compared with the 
particle period f2 -1 . In particular, the untrapped parti- 
cles can diffuse in J- variable and cross the separatrix 
k = 1 becoming trapped and vice versa. This diffusion 
may cover the phase space globally in contrast to stochas- 
tic layers around separatrices that usually develop in the 
case of quasi-monochromatic perturbations [ j35| . Coming 
back to the subject of this paper one can expect a simi- 
lar process that results in particle exchange between the 
trapped and untrapped regions as shown schematically 
in Fig. | 

In a general sense the origin of the background turbu- 
lence is not important in our simplified model. It is quite 
clear that there are many factors in the shock neighbor- 
hood that can drive such a turbulence and, hence, de- 
stroy the invariants of the particle motion. Our critical 
assumption is, however, that the particle diffusion in J- 
variable associated with this turbulence is slow compared 
with the regular motion: D(J)/Sl(J) <C Jf , where D is 
the quasilinear diffusion coefficient. On the other hand 
we also assume that the background turbulence provides 
sufficient mixing of particles in the region < z < L, i.e., 
DkoL/£l 3> J|, so that we finally impose the following 
constraint on the diffusion coefficient 
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1 D 



< 1 



(23) 



in the preceding section one easily finds 



Under these circumstances the escape flux upstream be- 
comes virtually independent of D i.e. of the downstream 
background turbulence. To explain the point consider 
the phase space in Fig. ^ again. As seen in the wave 
frame the shock front z = z s is escaping to the left at the 
speed uq = U2 — Ca ~ 112, where Ca is the phase veloc- 
ity of the wave propagating backwards in the local fluid 
frame as one excited by the backstreaming particles in 
the upstream region and transmitted then downstream 
p6| . Therefore, the particles above the upper branch of 
the separatrix S + , (k — 1) cannot penetrate upstream. 
At least some particles below S~ can in principle cross 
the plane z — 0. To identify them we first consider a 
particle having k(J) < 1 below S~. Its velocity v z os- 
cillates in time and depending on the invariants J and 
v the absolute value of particle velocity can exceed at 
least instantaneously the value uq and, therefore, such a 
particle can potentially take over the shock front. For 
this to happen it must reach z = 0, and since the diffu- 
sion in J is slow, it must come from a far downstream 
region. Therefore, it should exceed the velocity uo not 
only instantaneously but also on average (over f2 _1 ) to 
be able to get to the shock front and then to cross it. 
These particles are shown as the dotted area in Fig. ||. 
Therefore, the area between the separatrix S~ and the 
dotted area (where v z < — ito) should be virtually empty 
in the vicinity of the shock z = 0+, since the shock is 
effectively escaping from these particles and refilling this 
area by the trapped particles (above S~) and escaping 
(below S c ) due to the diffusion across S~ and S c is slow 
due to inequality (p3|). The same arguments can be ap- 
plied to the trapped particles which have v z ~ e > — uq 
and, therefore, cannot reach the shock. An exception 
should probably be made for some particles in the left- 
most trapped region. Their return upstream may occur 
when S~ crosses the line v z = —uq (not shown in Fig. 

provided that they have a proper phase in the wave 
while crossing the shock. However, this would be rather 
a reflection off the shock front than the leakage from the 
downstream region considered here. We shall not take 
this possibility further into account here (see Ref. j37|] 
for a reflection dominated injection scenario). From the 
above reasons we can identify the particles escaping up- 
stream with those below S c in Fig. ||. If the amplitude 
of the background turbulence satisfies condition (p3|), the 
flux of the escaping particles can be obtained from the 
ergodic arguments regardless of any details of their inter- 
action with the background turbulence. One may think 
of a quasi-linear plateau that also does not depend on 
the form of the wave spectrum as of a simple analog to 
this situation. For the particles on the curve S c we have 



1 

2^ 



V z dl\) = — Uq 



(24) 



' jfcK(k) 



— e 



- 1 



E(fc)-K(fc) 
J k 2 XK(k) 



(25) 



Equations ( p4| , p5| ) define a critical k = fc*(ito), so that 
particles with k < (dotted area in Fig. ||) escape 
upstream, whereas the rest are convected downstream. 



IV. INVARIANT MEASURE OF 
ADIABATICALLY ESCAPING PARTICLES 

In the preceding section we have defined a boundary in 
particle phase space that divides the shocked downstream 
plasma into two parts, escaping upstream and convecting 
downstream. In a four-dimensional phase space (v, z) 
this boundary is a hyper-surface given by the equation 
(see Eqs. (MM) 



V Z (J, A) = -Uq 



(26) 



Using Eq. (|T^) and the action-angle formalism introduced 



The unperturbed motion takes place on a 2-torus labeled 
by the two arbitrary chosen independent integrals of mo- 
tion, which are equivalent to two action variables. Reduc- 
ing the motion to an effectively one-dimensional we have 
used the action J as an action of the one-dimensional 
motion, and the integral A accounted for the motion in 
cyclic variables. These invariants were useful for the task 
of the two preceding sections. However, our final goal is 
to calculate the escaping flux starting from the Rankinc- 
Hugoniot relations, i.e. from the parameters of the up- 
stream flow, which unfortunately yield only the width of 
the downstream distribution in v, not its form, whereas 
the pitch angle distribution is implied to be isotropic. 
We thus need to transform the results of the preceding 
section to the variables v, y,, where /j, is the cosine of a 
pitch angle in the wave frame. A reasonable starting as- 
sumption about the particle distribution is that far down- 
stream where the thermalization of the plasma (also due 
to the interaction with the Alfven wave) is completed, 
the distribution becomes isotropic in pitch angle and, for 
example, a Maxwellian in v. 

As it was argued in the preceding section the back- 
ground turbulence will generally destroy the 2-tori and 
give rise to a relatively slow diffusion in A and J. If we 
assume in addition, that this turbulence is mainly due 
to the weakly dispersive Alfven waves or magnetosonic 
waves propagating almost (anti)parallel to the unper- 
turbed magnetic field and all in the same direction, then 
this diffusion is essentially a diffusion in pitch angle. In 
other words v remains invariant and the second order 
Fermi acceleration is not important. 

Our next assumption concerns mixing properties of 
the particle dynamics influenced by the background tur- 
bulence. In particular, we assume that the relevant 
phase flux is ergodic, and that the typical length scale 
L M of the pitch angle scattering satisfies the condition 
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1/fco <C <C L (see Eq.(p3|)). The ergodicity implies 
that when a particle wanders in the region z S (0,L), 
the time spent by it in the "escape" position is propor- 
tional to the size of the "escape" region. 

Thus, to calculate the fraction of the backstreaming 
particles that can cross the shock as a function of v 
we need to calculate the fraction of the hyper-surface 
v = const downstream from the shock occupied by the 
escaping particles. According to the preceding section 
this fraction can be represented as an invariant measure 
of these particles on the isoenergetic surface in the phase 
space as follows 



1 



(5 (v — v) dv'dfid(j)dz. 



(27) 



We have here introduced the spherical coordinates in the 
velocity space 



VfJ, 



vy/l — /i 2 cos 4> 
v yl — /i 2 sin <j> 



(28) 



The integration region T contains one wave period in z, 
i.e. 2tt and is confined in other three variables by the 
hyper-surface k — fc*(A,wo) = const to be obtained from 
Eqs. fl24| j25|) . According to the normalization used in 
Eq. ( |27| ) the full measure v(v) = 1, and Eq. ( p7j ) thus 
yields the fraction of particles escaping from the surface 
v = const. The boundary fc* of the region T is not 
a coordinate surface in the variables used in Eq. j27| ) 
over which the far downstream distribution is assumed 
to be uniform as discussed earlier in this section. To 
evaluate the integral (^), we therefore transform it to 
the variables already introduced in Sec. || as follows: 
(v, fi, cf), z) -> (A, q, r\, z), with 



2vy/l — fi 2 cos a + 1 + 2ev[i 




rj = VfJ, 



1 + v 2 (1 — p?) — 2^J\ — cos a 



(29) 



and a = <f> + z. The absolute value of the Jacobian of 
this transformation can be conveniently expressed after 
some algebra through the invariants of the unperturbed 
motion v and A 



d(X,r],q) 



d(v,/j,,<p) 
Thus the invariant measure 
1 



Av 2 TT 



J rewrites as 

5 (i>'(A, r), q) — v) XdX di] dq 



(30) 



(31) 



where v' 2 {X,rj,q) — 2H + (A — l) 2 , and H is defined by 
Eq. (|l5|). The next transformation which we perform is 



a canonical one, also introduced in Sec. |l[ i.e. (77, q) 
(J,ip)- Then we get: 



—v , S(v' — v) XdXdJ 
2v 2 J r 



(32) 



Transforming this integral to the k- variable introduced 
in Sec. |n[ after some simple algebra we finally obtain 



2V2 f k ' fa* 1 + y/ifc^-P + l 



.3/2 



k 3 J \k 4 v 2 - k 2 + I 



x K(k)dk 



(33) 



An escape of these particles is possible only for v > 1 and 
the escape probability behaves as is csc ~ kf (v — l) 3 ^ 2 for 
V — 1 <C 1, where k* is also rather small numerically for 
v < fth a j where w t h a m ay be characterized as a thresh- 
old velocity of the adiabatic escape and can be estimated 
from Eqs.(p4p4) as «th a ~ 1 + £■ For larger v, fe* rises 
sharply to reach fc* ~ 1. For v — > 00, v csc ~ 5 (l ~ 
The former formula for v csc reflects the shrinkage of the 
phase space of adiabatically escaping particles as v —> 1+ 
and does not mean that there are no other escaping parti- 
cles with v ^ 1 (see the next section) . The interpretation 
of the last formula is straightforward: these particles are 
not influenced by the wave and nearly a half of them 
escape. 

To obtain the distribution of escaping particles one has 
to multiply v csc by the thermal downstream distribution. 
To be specific we assume that the latter is a Maxwellian 
with the downstream thermal velocity vi- Then, using 
our dimensionless velocity in the wave frame v which is 
virtually the downstream velocity (Ca *C U2) the pith- 
angle averaged distribution of escaping particles can be 
written as 



2 

F esc (v) = r-fesc(v)/M(w)) 

1 - Uo/V 



V > 1 



whereas 



722 



TM 



■ exp 



2v 2 



(34) 



(35) 



with V2 = efcoVr/wci, and Vt = WT2/M is the down- 
stream thermal velocity. The factor (1—uq/v)/2 accounts 
for the limited fraction of the phase space at given v in 
which particles escape into the upstream half-space. Ac- 
cording to Eqs. (21 2^) we can estimate v 2 as 



v 2 < e— w V%e, 

U2 



(36) 



where the last value is valid for a strong shock with a 
compression ratio of 4. 

Since v csc starts to grow from zero only at v > 1 we 
infer that the contribution of the particles interacting 
with the wave adiabatically is exponentially small. 
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V. RESONANT ESCAPE 



it 



k u 2 



(42) 



Let us turn now to the resonant particles, that cannot 
be described by the simple formalism developed in Sec. 

since A 2 can be very small or even negative in this 
case. It is convenient to use the variables (v, /i, a) again 
(see Eq. (p8h), where a = <f) + z, and a new Hamiltonian 







H = yl — i-J. 2 cos a + -v[i 2 



which is related to A 2 through H = (v 2 
The exact equations then take the form 



(37) 
A 2 ) /2v. 



dn r, j . 

at 
da 
~dt 



&H 

da 



[i cos a 



M 



VjJL 



(38) 



The particle trajectories are shown in Fig. ^ as a contour- 
plot of the function A 2 with fixed v. Next we concentrate 
on the region of small and negative A 2 . First, consider 
the vicinity of the elliptic point of the Hamiltonian (|37j ) 
a = (mod 2ir), /i = /io. The equation for /zq reads 



M 



(39) 



and in the case of small e <C 1 and |/io| <C 1, the last 
equation reduces to a cubic one which yields for /io: 



Mo = < 



-2 |f| 2 sinh 
— 2f i cosh 
2f ^ sin 



i sinh 1 



i cosh 1 -%■ 

6 J 2 



isin" 1 



2tt 
3 



, < f < e§ 

, f > £3 



(40) 



where f = — 1). The bottom expression is strictly 
valid when f <C 1, otherwise a more accurate treatment 
of Eq. ( j39|) is needed. At the same time, for larger f 
the main contribution to the particle escape comes from 
the adiabatic region which was already considered in Sec. 
IV. Moreover, the downstream thermal particle distribu- 
tion falls off very rapidly in v, and the behavior of the 
escape probability at the lower energies is generally more 
important. We start, therefore, our consideration from 
the case v < 1. A particle that moves at the critical point 
a = 0, )jl = jjiQ can escape only when 



Mo v < 



u 



(41) 



(see Eq. (Pq)). Since |^o| is always small for £ < 
lfi ~ - ( 2s ) 1/3 for |f | < e 2 / 3 , and fi ~ 2e/3£ for 
|f | > £ 2 / 3 ), the last inequality can not be fulfilled for all 
resonant particles with v < 1 and therefore a threshold 
velocity D t h occurs. As it was argued in Sec. 
dimensionless velocity uo can be estimated as 



III, the 



where u 2 is the bulk velocity in the downstream region 
(see also Eq. ( |2l| ) ). In practice uq can deviate from the 
value ( (42| ) due to a number of reasons, e.g., due to a 
finite propagation speed of the wave and/or due to the 
fact that fco depends on vth- Therefore, we represent uq 
as uo = £ Ci where f ~ 1. Then, using Eq. ( |39| ) 
inequality (JO) rewrites as 



v > vth = CV (C + 1) + 



1 2 

2 +£ 



(43) 



To obtain the number of particles that are trapped into 
nonlinear resonance around the point a = 0, /i = /io 
and escape upstream we can use the same arguments as 
in Sec. IV in our calculation of the escape flux of the 
untrapped particles. Thus, given v > Wth we calculate 
the critical orbit around a = 0, (J, = (Xq for which 
the averaged velocity v z = vp, — — uq. Due to the anhar- 
monicity of the oscillations the averaged velocity v z that 
starts from the value v z = v/Iq at the bottom of poten- 
tial well, increases with the radius of the orbit, reaching 
finally — uq. Near the threshold where v — t>th <C 1, the 
condition v z = —uq is fulfilled for an orbit which is close 
to the critical point of the Hamiltonian. It is therefore 
convenient to introduce a new variable v = fi — fio and 
expand the Hamiltonian ( p7| ) at a = v = 0, retaining 
only cubic anharmonicity in \v\ <C 1, and neglecting a 4 
and a 2 vfiQ compared with a 2 « 1. Thus, the truncated 
Hamiltonian takes the form 



Hi 



a 

T 



2\ 2 MO 3 

Mo)^ ~Y" 



(44) 



We introduce again the action-angle variables (v, a) \— * 
(if), J), considering a as a momentum. The transforma- 
tion is generated by S = J adv. Thus, for ip, J we obtain 



J = 



ndv, 



da 



(45) 



where 



a = \J-no (v 3 - a 2 v 2 + a ) 

3 (i 

«2 = Mo 

2 VMo 



(46) 



and ao plays the role of the "energy" constant of the 
anharmonic oscillator (Q). It is worthwhile to perform 
the following transformation 



1 , 
v= ~a 2 (2z + l) 

and to rewrite a in Eq. ( f45|) as follows 

3 

2 



(47) 



d = | -a 2 j \/-Mo {z - z-i) (zq - z) (zi - z), (48) 
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where 



sin [ — sin 



2nn 



(49) 



Here a again denotes the "energy" of the oscillator ( |44| ) 
and varies between -1 and 1 in the potential well. Substi- 
tuting Eqs. (0,|||) into Eq. (|45|), we obtain after some 
simple algebra 



J = 3* 



2 (fc' 2 + fc 4 ) E(A) - fc' 2 (1 + k 12 ) K(fc) 
(fc' 2 + fc 4 )^ 



(50) 



where k 2 = (z - z_i) / (zi — Z-i) and fc' 2 = 1 - fc 2 . As 
it was mentioned, there exists a critical J = J , such that 
the particles with J < J c escape, whereas the particles 
with J > J e do not. To calculate J c we introduce the 
averaged p, as 



1 f 271 
fi(J) = — J dipfi(J, tp) 

Then J c will be defined by 

JX{J e )v = -uo = —Qe 
The last equation can be rewritten as 
v (P + /i ) = -eC 



(51) 



(52) 



(53) 



where v = (1/T) <f> vdv/aiv), a{y) is given by Eq. (46) 
and T = § dv / a. After a short calculation we obtain 



«2 



Vfc' 2 + fc 4 



~ (2 - fc 2 + \Jk 12 + fc 4 ) 



E(fc) 



K(fc) 



(54) 



Inserting the last expression into Eq. ([53|) we first obtain 
fc„ as a root of the equation 



v {v(K) + fi ) = -e(, 



(55) 



and using then Eq. (|50|), we finally obtain J e = J(fc„). 

Normalizing the full measure of the particles at v 
const to unity 



1 

47T J_ 1 



dft 



da = 1 



(56) 



and adopting arguments similar to those already used in 
Sec. IV, we find for the invariant measure of the reso- 



nantly escaping particles 



^esc . 

47T 



dadu 



J<J B 



Transforming dadfi — > dipdj we thus find 



1 



1 



J(fc*) 



(57) 



(58) 



Not far from the threshold velocity (v > fth), where fc* 
is small, from Eq. (55) we obtain 



16 



(I2K 



-Mo 



V 



(59) 



Since 



J(fc) 



3* (^-e) 



5/2 



2* 



Mo 



for 0SC we have 
V2 



n csc « ^ (3^ 2 - 2(« - 1)) 3/2 ( 1 + — 

which in the case 1 — v > e 2 l z simplifies to 
2 



3v 



(l + C)(l-u) 5 (u-vth) 



(60) 



(61) 



As v grows approaching unity, Eq . (p9 ) becomes invalid 
and in the case opposite to Eq. (|59|) i.e., when fc' -C 1 
from Eq. (pH), we find 



1 — 16 exp 



-3v- 



2(v - 1) - 3 M 2 



2(u- l)u + 3£/i C 
1 - 16cxp|-3 • 2 1/3 e- 2/3 /c} 



(62) 
(63) 



Note, that the last expression is valid only for |1 — v\ < 
e 4 / 3 . The escape probability then takes the form, Eqs. 



5V2tt 



5/2 



Mo 



(64) 



The dependence £l csc (v) as calculated for small fc and for 
fc ~ 1 is shown in Fig. [|. Since the downstream thermal 
distribution falls off very rapidly the resulting escape 
spectrum will have a maximum in v close to the point 
v = Vth ~ \- The distribution of the escaping particles 
can be written again as 



1-Uq/v 



^esc/M(l>) 



(65) 



For larger v ^ 1, when the approximation \/j,q\ <C 1 
breaks down, formula (|64|) becomes invalid as well, and 
a more accurate consideration of Eqs. ( |38|) is needed in 
this case. At the same time f2 0SC ~ e for v ~ 1, as 
may be seen from Eq. (|64|), and the escape in the region 
f — 1 > 1 is dominated by particles that interact with the 







wave adiabatically. Thus, for v > 1 the distribution of 
escaping particles can be given by Eq. (|34|) . It should be 
born in mind, however, that the maximum of the escape 
distribution is rather close to v = i for small e, so that 
formulae (|6l]) and (^i|) provide virtually the kernel of the 
distribution of escaping particles, whereas Eqs. (HHJ) 
describe the tail of this distribution. 



them at the current level of description. The simplest as- 
sumption is that upon crossing the shock these particles 
behave more or less like the protons. In other words, their 
excess (over u 2 ) velocity, i.e. u\ — u 2l is spread around u 2 
and we assume that Vr a — Vt- For the purpose of sim- 
plicity and for extracting the dependence upon A/Z we 
also assume that the thermal distributions of all species 
a is equivalent to that of the protons 



VI. INJECTION EFFICIENCY VERSUS MASS 
TO CHARGE RATIO 



fa 



n 2a 



(27T) 



3 / 2 „,3 



■ exp 



2a 



2?j 2 

LV 1a 



(67) 



One important aspect of any injection mechanism 
should be its dependence upon the mass to charge ra- 
tio of different species. It is obviously so in the mecha- 
nism suggested in this paper. Indeed, in the case e <C 1, 
the leakage upstream must be controlled by the param- 
eter fcop Q , where p a = (Vr a / ijJ ±)(A/Z), is the Larmor 
radius of a species a, i.e. Vr a is a corresponding ther- 
mal velocity, A and Z are the mass and charge numbers, 
respectively. It is clear that strongly magnetized parti- 
cles (kop a *C 1) cannot be injected, whereas unmagne- 
tized particles (kop a ^> 1 are injected as readily as in the 
case without magnetic field. For protons this parame- 



ter is fcoPp 



t'2 



> 



e almost by definition, simply due to 



the fact that both the wave and the thermal distribution 
downstream originate from the same upstream flow (see 
Eq. (|36])). To confine particles effectively the parameter 
e must be rather small but not too small- otherwise the 
upstream turbulence cannot be exited by a weak proton 
beam. This means that protons are close to a watershed 
between the species that cannot be injected by this mech- 
anism (these are apparently electrons, only) and particles 
with higher A/Z whose injection efficiency increases. 

According to the preceding section the most important 
physical quantity that regulates the escape flux is the 
threshold velocity Vth which is the same for all species, 
Eq. (ff3|), provided that in the definition of the normal- 
ized velocity v, Eq. (||), one substitutes uj± a = uj±-(Z/A) 
instead of Thus, the subjects for injection are only 
the particles with v > Vth — 1/2, or (kov/oj±)(A/Z) > 
1/2 for unnormalized v. The escape probabilities v esc (v) 
and Q csc {v) (Sees. [V V]) as functions of dimensionless 
velocity are also the same for all sorts of particles. The 
quantity that discriminates particles against A/Z ratio 
in dimensionless variables is obviously the downstream 
thermal velocity V2, Eq. (p6[). For a species a we thus 
have 



Then, we may calculate the density |38| of injected par- 
ticles 



2dv 



v.<-u 1 - u a/v 



fa^c 



(68) 



where f2 csc = Q csc for v < 1 and fl csc — v csc for v > 1 (see 
Eqs. ( |34| ) and (^5|)). The last equation can be evaluated 
to 



n taj = \/~-§^/ v 2 dvn csc (v)exp 



J 2a Jv 



2< 



(69) 



where v = max {vth, uoa} uo a = uqA/Z. According to 
Eqs. ( |60|j6^ ) and Fig. [|tne function Cl esc (v) rises sharply 
from r2 e sc(vth) = to become approximately constant 



o 3 ^ 
U = — — e 

07T 



in the region t>th + ^ v ^ 1, where a = 3 5 / 2 /107r. For 
v > 1 the escape flux is dominated by adiabatic particles. 
To simplify the algebra we substitute Cl csc ~ fl into Eq. 
( |69| ) and shift the lower limit w t h — ► w t *h = w *h + where 
a ~ a and extend the integral to oo. Starting from v = 1 
we may use the high energy asymptotic result v csc = (1 — 
l/v)/2, multiplied by 1 — 2Hq to compensate the above 
extension of the contribution of resonant particles. This 
simple interpolation yields for the protons, rj = «f n j/«2p 



( u t* h 



VV2 



('2 



2 ft 

7T V 2 



(70) 



where 



l'2n 



v 2 - 



Vra A 
V T Z 



(66) 



$(x) = —= \ e" £ dt 

V 7 !" JO 



We may use now Eqs. (|3J) and ( p5| ) for calculating the 
distribution of escaping particles of a sort a upon sub- 
stituting V2a instead of v 2 - It is of course assumed that 
all the arguments of Sec. Ill are valid for these particles 
as well. There exists, however, the problem of the ther- 
mal velocities Vr a ■ It is indeed very difficult to quantify 



For other species one may obtain a similar formula from 
Eq. (|69|). We illustrate it by plotting the injection ef- 
ficiency r\ a normalized to the proton efficiency. This is 
shown in Fig. for v 2a = v 2 A/Z and a = 1.5c. 
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8tt 



A 



-n b mpCAWb 



(71) 



We have merely introduced an additional factor A ~ 
Avu/vb < 1 where Auii is the beam width in parallel 
velocity. This factor appears because the beam relax- 
ation occurs under the constraint of conservation of the 
(zero) particle flux J v\\fbdv\\ on a given diffusion line 
Vj_ + (vh +Ca) 2 = const in velocity space rather than un- 
der the conservation of the phase density f fbdvu along 
this line. The reason for such a factor A may be un- 
derstood from the observation that particles starting to 
escape at velocities ~ — Aun, while being turned around 
can hardly acquire positive velocities that are appreciably 
larger than +A«ii before returning to the shock. Again, 
because the above mentioned particle flux must be zero. 



VII. A SIMPLIFIED SELFCONSISTENT MODEL 

So far, we have considered particle escape under a pre- 
scribed wave spectrum downstream. However, as we em- 
phasized, this escape mechanism ought to possess a very 
distinct selfregulation. Indeed, there is a strong negative 
feedback between the wave intensity and the density of 
the escaping beam- if the beam is weak and excites thus 
only weak waves, the leakage will be increased to produce 
stronger waves. Similar arguments lead to decreasing 
the leakage if the beam is too strong. Therefore, both 
the beam intensity and the turbulence amplitude must 
rest at some definite and unique level. What makes this 
situation differ from the standard quasilinear theory of 
beam relaxation in homogeneous plasmas is that "beam 
relaxation" here means actually its return to the shock 
front via the cyclotron interaction with the self-excited 
MHD waves (see Ref. [[l3| for a detailed description of 
this process). Hence, a plateau does not form in fact 
and the relaxation length means simply a distance at 
which the majority of beam particles are turned around 
and swept back to the shock. The nonlinear wave phe- 
nomena are assumed to be unimportant, which implies 
that the corresponding time scale tnl > Ir/ui. Thus, 
the fraction of the beam energy that may in principle be 
channeled into the plasma heating through the nonlinear 
wave-particle interactions is correspondingly small and 
the wave energy at the shock front may be calculated 
from a simple energetic balance. 

First we note that the beam energy is mnbV^/2, where 
I ft> | ~ u i is the beam velocity in the upstream frame and where 
rib is its density. Since it scatters back quasi-elastically, 
namely around scattering centers that move at the low 
velocity —Ca, only a Cf^/vb *4C 1 fraction of beam en- 
ergy may be converted into waves. A complete quasilin- 
ear theory of cyclotron beam relaxation in homogeneous 
plasmas has been developed in Ref. p9| . One can show 
that the expression for the wave energy released by an 
unstable beam as calculated by these authors is also ap- 
plicable for the case considered here. Thus, for the mag- 
netic field perturbation upstream we may write 



Since the downstream field B± — rB± u , where r — u\lu-> 
is the shock compression ratio we may rewrite Eq. (|7l| ) 
as follows 



1 



1 



-r kMpj] 



(72) 



where Ma — Vb/CA and 77 = nb/772 is given by Eq. (|7C 
Since 77(e) is a monotonically increasing function, Eq. 
( f72| ) determines a unique value of e and thus a unique 
injection rate 77. 

Consider first the case Ma 3> 1. To the leading ap- 
proximation in e <C 1 and substituting V2 = Voe, r = 4 
into Eq. (170) we obtain 



6 



77 = 



^T^exp 



V* 2 

v th 

"6e 2 



Eq. (|72|) then rewrites 



1 



_=CM A exp(-^ 



(73) 



(74) 



where 



C 



192 



v : h A 



57r 3/2^ 

For CM a 1, and assuming C ~ 1 we thus obtain 



L A = In M A 



(75) 



y th 



"6 In M A 



For the injection rate 77 we have 

3£a 



V 



16AM A < h 2 



(76) 



One may put here v* b ~ 1/2. It is seen that the injection 
rate formally vanishes with I /Ma as 77 ~ M^ In Ma- 
However, already for e < 1/4, B± u ~ B z and the esti- 
mate of the wave amplitude might need some correction. 
Clearly, the above 77 scaling is not applicable in the limit 
B z — > in which 77 vanishes. For such a weak magnetic 
fields different mechanisms of beam relaxation must be 
considered. 

For moderate values of e < 1, Eq. ( |72] ) can in principle 
be easily solved numerically with 77(e) given by Eq. (|7Q|). 
We know already that the above formulae are applicable 
for rather small values of e < 0.2. On the other hand 
e cannot be too large in any case, in fact it cannot be 
larger than about 0.4 to satisfy Eq. (f72|). Thus, what is 
actually needed is a reasonable but simple approximation 
of 77(e) in Eq. ( |72] ) in the interval 0.2 < e < 0.4 to resolve 
Eq. ( [72] ) for e. It is convenient to use the following 
approximation 
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1 - C 2 £ + C 3 E* 



which is shown in Fig. || for c\ = 0.0105; c 2 = 4.91; C3 = 
6.58. Denoting q = 32AA/aci we obtain the following 
solution for e as a function of q 



C'2 



2(c 3 - q) 



1 - 



K c 3 -9) 



(78) 



It is shown in Fig. fj] together with 77(g). One sees that 
the injection rate depends rather slowly on Ma as in the 
case of higher Ma (smaller e) considered earlier. 



VIII. THE LEAKAGE PROCESS IN A NUTSHELL 

The mechanism of ion leakage considered in the pre- 
vious sections unfortunately requires more calculations 
than seems to be appropriate to its physical simplicity. 
At a phenomenological level this mechanism is almost as 
simple as the escape from an oblique shock where for a 
particle to catch the shock it must move at the speed 
U2/cos6„_b 2 — u 2/e along the field line. Here 6„s 2 is 
the angle between the downstream magnetic field and the 
shock normal). Such an escape has been extensively stud- 
ied in Ref . |39|] . Even if the shock is quasiparallel but the 
magnetic field is locally oblique to the shock normal most 
of the time, the same kinematic escape condition holds for 

most 



magnetized particles. As we have seen in Sec. VI 



of the protons must be magnetized and in the case of the 
turbulence dominated by a circularly polarized Alfven 
wave the lowest energy particles that can escape have a 
velocity v — Vth — 1/2 («th — u 2 /2e in unnormalized 
variables). According to the phase plane shown in Fig. || 
these are the particles that move towards the shock being 
close to the point where a = <fi + z = 0; fi = (1q < 0. 
Therefore, they fall into the cyclotron resonance with 
the wave (4> + z — 0, <fi — —z = —/iQvt but they spi- 
ral as electrons in the unperturbed field B z trying to 
follow the magnetic field line and to minimize thus the 
Lorenz force, the only force in our model that can pre- 
vent their escape. The inclination of their orbit to the 
z~ axis is, however, about a half of that of the magnetic 

field: v±/v z = ^Jv% + jv z ~ (1 — v)/e ~ l/2e whereas 

B ± /B z = 1/e. 



As we have seen in Sec. VI, in the case of very small 
values of e these particles make the bulk of the leakage. 
Since they are concentrated in a relatively small region 
of the downstream phase space, this allows us to fore- 
cast their energy and angular distribution, just as they 
appear upstream. First, their energy per mass in the 

A 



downstream frame must be somewhat above E = 
where «th — 1/2 in the dimensionless variables. Fur- 
thermore, since v±/v z ~ l/2e ^> 1 the energy of leaking 
particles is predominantly in the perpendicular motion. 
In the unnormalized variables we thus have 



2 B " th " 



whereas the parallel energy may be estimated as £» ~ 
Ae 2 £^ . It should be noted that the pitch angle scattering 
upstream may change this relation to a certain extent. 
On the other hand it is in a reasonable agreement with 
the results of the strong shock simulation |) 23 1 . 



IX. INJECTION EFFICIENCY. COMPARISON 
WITH HYBRID SIMULATIONS 

As we emphasized in the Introduction the calculation 
of the flux of leaking particles alone does not solve the 
problem of injection. The main result of injection the- 
ory should be the high energy asymptotics of a spectrum 
that emerges in a steady state when the leaking particles 
repeatedly cross the shock and achieve energies sufficient 
for describing them by the means of the standard theory 
of diffusive shock acceleration (see e.g., ref. |lq|). The 
mathematical formalism of injection theory has been de- 
veloped in Ref. |13|. Now we may apply it to the dis- 
tribution of leaking particles (thermostat distribution) 
calculated in the present paper. We also have to bear in 
mind that particles that cross the shock more than once 
(higher generation of injected particles, beam 2, etc., see 
Fig.|l|) are still subject to the filtering on their way back 
upstream due to the interaction with the downstream 
trailing wave. As we have seen this interaction weakens 
with the energy and the thermostat becomes transparent 
to particles with v 3> ujj_/ko (unmagnetized particles). 
With this in mind, the whole algorithm may be outlined 
as follows. 

Suppose that some fraction of the downstream plasma 
leaks upstream to form at z — the one-sided distri- 
bution F(v),«|| = v' z < (see Fig. |lj). Here v' is the 
velocity in the shock frame and we keep our notation v 
for the wave frame in the downstream medium (almost 
the downstream frame). Due to pitch angle scattering 
in the upstream medium these particles turn around and 
eventually cross the shock in the downstream direction 
forming the distribution F + (y), v' z > which can be 
written as F + = L\F, again at z = 0. The linear oper- 
ator Li, the upstream propagator, can be obtained from 
the solution of the kinetic equation fl3||. According to 
the thermostat model in use (Sees, p A||IIID these parti- 
cles penetrate further downstream through the thermo- 
stat mixing up with the hot downstream plasma. At the 
same time they are pitch-angle scattered on the back- 
ground turbulence so that some of them acquire negative 
velocities and move back to the shock. We denote their 
distribution within the thermostat by F~ . For F~ we 
thus have 



F~ = L 2 L X F + /] 



M 



(80) 



Here L 2 is the downstream propagator, and /m is the 
distribution function of the downstream thermal plasma 
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that emerges upon the first crossing of the shock interface 
(without higher generations). In Sees. |v|.|v]we assumed 
for simplicity that /m is a Maxwellian distribution so that 
L 2 /m ~ /m because it is isotropic in the wave frame. 
Now the calculation of the injection spectrum that ap- 
pears just upstream of the shock is nothing more than 
the calculation of the spectrum of leaking particles al- 
ready made in Sees. I V,V] with Jm in Eq. ( |65| ) replaced 
by F~ from Eq. (|S0|). Thus, the distribution of injected 
particles for F takes the form 



F = tL 2 L x F + r/i 



M • 



Here the function t(v) is given by (see Eq. (pq) ) 

2 



t{v) 



1 - U /V 



(81) 



(82) 



and may be interpreted as a thermostat transparency co- 
if there 



1 as v — ► oo. 



efficient. According to Sec. IV - 
were no scattering upstream (L\ = 0), Eq.(|3lj) would 
be equivalent to Eq. (^5|) that yields the distribution of 
escaping particles given the thermal distribution down- 
stream, /m- In general, Eq. ( |Sl| ) is an integral equation 
for F. The kernel L 2 L\ has been calculated in |l3| where 
also the solutions of this equation have been studied for 
t = 1 and various functions /m which should have mim- 
icked the effect of thermostat filtering (r < 1). 

Expression (82) implies only an adiabatic leakage 
which is appropriate for particles with v > 1, Sec. IV. 



If the resonant particles leak as well, the escape proba- 
bility f2 esc from Eq. ( |65| ) should be added to v csc in Eq. 
( |S2| ) and it will dominate in the region v < 1. However, 
unlike the leakage of adiabatically interacting particles, 
the leakage of the resonant particles is very sensitive to 
the wave polarization. This may be understood from in- 
spection of Fig. H drawn for an A-wave. The MS-wave 
case can be obtained by flipping the phase portrait in 
Fig. m since the MS polarization corresponds to e < 
and the Hamiltonian ( |37| ) is invariant to the transfor- 
mation e — > —e, (i — > —fx. Thus, the candidates for the 
resonant leakage from an MS-wave would be the particles 
marked by 2 and 3 and them alike, i.e., those circulat- 
ing around a fixed pont at a = —n (mod 27r), /i > 
in Fig. ||. However, they have relatively low values of /2 
and rough estimates show that they cannot escape. At 
the same time, according to our discussion of the leakage 
from the thermostat in Sec. Ill such particles can po- 



tentially escape from the region immediately behind the 
shock front provided that the wave field is sufficiently 
perturbed. This point should be born in mind when we 
compare our results with hybrid simulations below. 

Most of the hybrid simulations are essentially time de- 
pendent, since the shock runs through a finite spatial do- 
main. Eq. ([y]) implies a steady state and for this rather 
preliminary comparison we select only a hybrid simula- 
tion [p3[ where the simulation box was anchored on the 
shock front and a quasi-stationary spectrum was devel- 
oped. A more thorough comparison with other numerical 



results will be done elsewhere. As we have seen, the most 
important parameters that determine the distribution of 
leaking particles are the amplitude, the wave number and 
the polarization of the trailing wave. The injection spec- 
trum is then formed depending primarily on the shock 
compression and to some extent on the spectrum of the 
background turbulence upstream and downstream that 
enters the propagators L\ and L 2 in Eq. ( jsi"| ) ijl^JfOf . 
For the purpose of comparison we assume the MS po- 
larization as observed in |^3| (see, however ref. [^J) and 
therefore discard the contribution of resonantly escaping 
particles. 

All other quantities needed for calculation of the in- 
jection spectrum can be obtained from the above results 
given Ma and Ms and from RH conditions. However, 
shocks that form in simulations do not follow the latter 
exactly, due to the losses through the boundaries of sim- 
ulation box and other reasons discussed earlier. There- 
fore, it is appropriate to take some critical parameters 
directly from simulations. For example, the total com- 
pression ratio obtained in |^3| is close to 4.2 exceeding 
the limiting value of 4. At the same time the shock is 
noticeably modified, so that the local compression ratio 
is substantially smaller. We make our comparison taking 
r = 4.0 which also corresponds to a strong shock. The 
downstream temperature slightly deviates from the RH 
prescriptions as well and we take it from the simulations 
(T = 2 • 10 6 K) in order to ensure coincidence of the ther- 
mal parts of the spectra. The amplitude parameter e for 
the shock of Ma = 5.25 in |23| may be calculated using 
Eq. (|72|) or @. We estimate A = Auy/i; b w 1 which 
is appropriate for a strong shock, and obtain e ~ 0.3. 
This e is in reasonable agreement with the simulation re- 
sults. It should be reminded however, that Eq. (|7^) was 
derived for the resonant leakage and a somewhat differ- 
ent although very similar equation should have been used 
in the case of adiabatic leakage which would have given 
a somewhat higher e. However, some of the resonant 
particles probably leak in the simulations and we adopt 
Eq. (78) for our estimate of e. Furth ermo re, in the self- 
consistent determination of e in Sec. VII only the "first 
generation" particles are taken into account which clearly 
leads to an overestimated e in Eq. (|72]). 

We calculate the wave number ko from the condition 
of frequency conservation across the shock transition, 
ko = k u r(MA — 1)/(Ma — \/r), where fc u is an upstream 
wave number that we in turn obtain from the resonant 
condition k u Vb ~ w c i that we wrote here to the leading 
order in I/Ma- In general, the diapason of unstable wave 
numbers may be quite broad since the escaping beam is 
broad in vu . An accurate calculation of the most unstable 
A: is a difficult problem since the beam distribution de- 
pends on this k as well and we restrict ourselves to a sim- 
ple estimate based on the mean beam velocity. Namely, 
the escaping beam occupies in velocity space at least an 
interval —Vth a^j_/^o + u 2 < V\\ < where v» is a dimen- 
sional particle velocity in the shock frame and Vth a is a 
demensionless threshold velocity for an adiabatic escape 
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(see text below Eq. (82)), v t h a ~ 1 + e roughly indepen- 
dent of the wave polarization. From this considerations, 
we obtain u>±/koUi all 

We compare our analytic calculations (see Appendix 
for more detailes) with the hybrid simulations p3| in 
Fig. [|. The downstream Maxwellian is the same in both 
cases and it is drawn with the thin line. The squares are 
from simulation whereas the heavy line shows the result 
of integration of Eq. (|Sl|, Al). The slope of the energy 
spectrum oc E~ a , with a = 3/(r— 1)«1 that must form 
at energies sufficiently higher than the thermal energy 
according to the standard theory of diffusive shock ac- 
celeration is also shown by the dotted-dashed line drawn 
at arbitrary height. Finally, the dashed line shows the 
solution of Eq. (81) with r = 1, i.e., for a completely 
transparent thermostat pi] . 

The first conclusion that may be drawn from Fig. |^ is 
that the effect of particle filtering by the thermostat is 
very strong and reduces the injection rate by one order 
of magnitude compared to the case of "free" injection, 
i.e., without the strong wave particle interaction down- 
stream. Furthermore, the agreement with the simula- 
tion spectrum is very good, although the latter does not 
exhibit a correct high energy asymptotics, most proba- 
bly due to the losses of high energy particles. It is rea- 
sonable to assume that if there where no such losses a 
correct spectral slope (with probably somewhat higher 
amplitude) would be achieved in the simulation at an 
energy where the slope has now a minimum (intersection 
point with the analytical spectrum) . This means that the 
two methods would produce similar injection rates if the 
latter is understood as an amplitude of the high energy 
asymptotics oc E~ a . The lack of low energy particles in 
the analytical spectrum should be attributed to the ab- 
sence of resonantly leaking (or reflected) particles that 
are probably still present in the simulations. This might 
also slightly underestimate the injection rate as it is per- 
haps the case in the simulations because of the losses. At 
the same time we feel that whatever physical ingredients 
(like the resonantly leaking particles) are added to the 
above calculation scheme it will not change the injection 
rate significantly. This is due to the strong selfregula- 
tion of the leakage (injection) process. The shock seems 
to leak always at a critical rate which is just enough to 
confine the downstream plasma through wave generation 
as it was discussed earlier. This fundamental aspect of 
particle injection at quasi-parallel collisionless shocks has 
been foreseen by previous authors, e.g., pj,p|,^2|. 



X. OTHER POSSIBLE APPROACHES, EXISTING 
AND PROSPECTIVE 

Existing theories of shock dissipation and shock accel- 
eration have not included the injection of the shocked 
plasma into the foreshock region selfconsistently. An 
important insight provide hybrid simulations, but being 



substantially limited in space, time and particle energy, 
they miss the backreaction of accelerated particles on the 
shock structure and, therefore, on injection and shock 
dissipation. Monte Carlo simulations ( e.g., (i^]) include 
the backreaction but they completely ignore the feedback 
from the turbulence excited by injected particles them- 
selves which may reduce the injection rate by an order of 
magnitude without major changes in the flow structure. 

Although the collisionless shock phenomenon is very 
complicated, the necessary information about the source 
of leaking ions can be inferred from the ordinary jump 
conditions. The latter show that at least behind a strong 
shock the plasma distribution is so broad that a very large 
fraction of it can in principle escape upstream. This will 
certainly smear out the subshock as a distinct structure, 
unless this thermal return is choked by fast unstable cou- 
pling with the incoming flow. Thus, the problem is not 
the source of the particles to be injected upstream, but 
rather the opposite, i.e. how to confine most of them on 
the downstream side of the shock and, of course, how to 
calculate the distribution of the rest which is leaking. 

One may attempt to do this in several ways. Firstly, 
one can invoke the processes occurring at the very shock 
interface and immediately behind the shock, before the 
thermalization of the plasma flow is completed. One 
obvious candidate is the electrostatic barrier appearing 
e.g., at the subshock in the high Mach number hybrid 
simulations || (some further discussion can be found 
e.g., in a review paper ]44[]). This is critically important 
in the quasiperpendicular shock mechanism [Iq| . Quest 
H noticed, however, that its impact on inflowing ions in 
strong quasiparallel shocks is very modest due to the re- 
markably perfect compensation of the electrostatic force 
with the Lorenz force. On the other hand, this conclusion 
may not be true for the backstreaming ions. 

The next possibility consists in the already mentioned 
subshock reduction by a pressure gradient built up by 
the intense return beam in front of the shock. This is 
the mechanism at work in Monte Carlo simulations jl3) . 
In fact, it is the only mechanism of selfregulation of the 
injection process at the subshock level in Monte Carlo 
models since they operate under a prescribed scattering 
law. It is, however, not to be confused with the process 
of a large scale shock modification by diffusively accel- 
erated high energy particles which can also reduce the 
subshock. Although these two processes of shock mod- 
ification are physically very similar, the latter operates 
over a much larger spatial scale and depends, besides the 
injection rate, on factors that have nothing to do with 
the subshock physics, like losses at highest energies [ ^6| . 
Note that the selfregulation mechanism suggested in this 
paper works very efficiently regardless of (or along with) 
the flow modification in the precursor and the above- 
mentioned feedback from the high-energy particles. 

The simplified model considered here gives explicit 
formulae for the distribution of backstreaming particles 
upon the wave amplitude (through the parameter e). 
The latter has, in turn, been calculated by considering 
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the transformation of beam energy into the wave energy. 
This upstream wave, driven by the unstable beam, may 
also be a subject to one of the known saturation mecha- 
nisms, especially in the case when a compressional com- 
ponent of the wave field is added. Besides the quasilin- 
ear beam relaxation considered here, these may be wave 
steepening p7| , other processes of nonlinear wave trans- 
formation, or beam trapping (see e.g., Ref. |29|). Nev- 
ertheless, like in the case considered here, the wave am- 
plitude should be calculated as a function of the beam 
intensity, providing the injection efficiency with no pa- 
rameterization. 



XI. LIMITATIONS 

We have assumed the unperturbed particle motion to 
be determined by a monochromatic wave which is an ex- 
treme idealization in the shock environment. At the same 
time it is often the case in wave-particle interaction, that 
the form of the wave field is not important for particle 
trapping- only the depth of a 'potential well' is impor- 
tant. The worrying situation here might be created by 
resonantly escaping particles (Sec. 0), that seem to re- 
quire a 'fine tuning' in wave-particle interaction. On the 
other hand, the background turbulence that was assumed 
to be sufficiently strong (Eq. (]23|)), certainly diminishes 
the role of this fine tuning by destroying the integrals of 
regular dynamics. Put another way, a particle escapes 
not because it stays in exact resonance with the wave 
for a long time which would hardly be possible for any 
realistic wave field at a shock, but because it appears at 
the right place in phase space while being close to the 
shock front. Otherwise its motion may be quite irregu- 
lar. The key element of our treatment that allowed us to 
calculate the escape flux under a restricted knowledge of 
chaotic particle dynamics was, of course, the ergodicity 
assumption. 

The injection scheme presented here will require cer- 
tain modification when applied to the case of finite # n B- 
Scholer et ai, demonstrated by means of hybrid simula- 
tions |?7j , that the contribution of particles staying suffi- 
ciently long at the shock front is increasingly important 
in this case. Within our scheme, these particles can be 
formally identified with the particles that are in nonlin- 
ear resonance with the trailing wave, have the averaged 
velocity v z w — uq (marginal escape) and are close to the 
shock front. Upon interaction with it they gain energy, 
although the mechanism wherby it hapens is yet to be 
studied. 

Generally, the tight link between the escape flux and 
the wave amplitude emphasized in this paper is the 
essence of the selfregulation of the shock dissipation pro- 
cess given e.g., by Eq. (|7l|). However, depending on the 
age and size of the shock this may be not the only way the 
shock regulates its own energy dissipation and particle 
acceleration. In nonlinearly accelerating shocks the sub- 



shock strength may be significantly reduced. Also, the 
deceleration of the flow in front of the shock by high en- 
ergy particles drives the subshock Mach number to lower 
values which may have an important impact on both the 
injection and the overall flow structure near the shock. 



XII. CONCLUSIONS AND DISCUSSION 

We have demonstrated that the large amplitude wave- 
train can efficiently filter the warm downstream plasma 
in its leakage upstream, scattering back typically no more 
than 5 percent of the downstream protons in the case of 
a strong shock and a left-hand polarized (Alfven) wave. 
MS-type of polarization results in noticeably better con- 
finement of the hot downstream plasma and weaker leak- 
age (injection). This means that the MS-turbulence is 
more suitable for maintaining a distinct quasi-parallel 
shock structure than the A-turbulence. The spectrum of 
high energy particles accelerated out of the backstream- 
ing beam is calculated with the help of injection theory 
p3[ . The resulting spectra are: (i) in reasonable agree- 
ment with the broad dynamical range hybrid simulations 
to date [p3| (ii) they evolve into a standard power-law at 
higher energies (iii) their intensity may easily exceed the 
threshold of the nonlinear acceleration regime (see be- 
low). 

It is needless to say that a reliable calculation of injec- 
tion rate, that take account of all essential interrelations 
between physical processes like the leakage/reflection, 
wave generation, particle trapping and shock modifica- 
tion by energized particles, could dramatically improve 
our understanding of how strong shocks accelerate par- 
ticles. Recent analytic solutions [[l6 48 for nonlinearly 
accelerating shocks (i.e. shocks whose structure may be 
almost entirely determined by accelerated particles) show 
that the dependence of the acceleration efficiency upon 
the injection rate has a critical character allowing for 
extremely different solutions at quite close or even the 
same injection rates. Therefore, the studies of the en- 
ergetic particle (cosmic ray) production ^] in such 
shocks or, in other words, of how the shock distributes 
its energy between thermal and nonthermal components 
of the shocked plasma should perhaps be focused on the 
subshock where particles are injected into the accelera- 
tion process. The necessary subshock parameters should, 
however, be determined selfconsistently from kinetic non- 
linear calculations of the shock structure like those men- 
tioned above. 
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APPENDIX: 

In an expanded form Eq.(j8l|) can be written as (see 
L3| for further details) 



3 1 /3r 2 (2/3) c 



T J ridrjAi k/ — r) 



Ai (k 1/3 ^v) ' F(vi) + Tf M (v) (Al) 



Here the following notations have been used 
C± = v± u 2 , 



= \/v 2 — 2 Au (C+v + ifrfj), Au = u\ 



U2, 



ip = \Ju\ + C+C- - 2A<+ 17 - u u 



rjo = min(l,C-/2Au), 

r is the Gamma function, Ai denotes the Airy function, 
and fj ~ 2/3. 

The coefficient K\ is a certain functional of the spectral 
density of the background turbulence downstream, that 
ensures pitch angle scattering and have been discussed 
already in Sec. III. We take m from JTa], Eq.(67), using 
slightly different notations 



where 



and 



2 cosh 



cosh 1 (a 



1) 



1 



(A2) 



r(i/3) 

3 4 /3r(2/3)A^ /3 



C-Ja 



C- 



dC(i - C/C-Y/d(-0 



Here D is the diffusion coefficient in velocity space nor- 
malized to its value at £ = as a function of the resonant 
(C ~ u 2 + oJd/k) parallel velocity of the ions in the shock 
frame. 

The coefficient K has a similar meaning as k\ but it is 
related to the particle transport in the upstream medium. 



k = —protect |3 — 7 + 2^/7(6 — 7) 
1 



sinh 



sinh 



7 2 + 9 7 - 27/2 \ 
V7(6~7) 3 / 2 J 



(A3) 



where 



and 



4^ 3 V3 



27r 6 (2/3)A 2 



A 



1 



dc(i - c/cr/D(-o 



Here — V — u\ is an upstream analog of calculated 
for the upstream absolute value of particle velocity V. 
Note that in |l3|] the expression for K\ (Appendix C) has 
been given erroneously only for the case 7 > 6 (even for 
7 ^> 1) and the formula for 7 contained a misprint. At 
the same time the case 7 < 6 has been actually considered 
in numerical examples, however, with the correct numer- 
ical values of n\ and 7. Generally, there still exists some 
arbitrariness in choosing the parameters k and k,\ since 
the details of the spectra of background turbulence are 
not determined in injection theory Jl3| ] . Nevertheless, the 
resulting particle spectrum may be calculated because it 
is rather insensitive to parameters k and k%, although 
they influence slightly the slope of the spectrum at high 
energies. In example given in Sec. IX we put k = 1.3, 
Ki = 1.15, ignoring their possible dependence on particle 
energy. Note, that the case k w Ki « 1 corresponds to 
a simple assumption D = const which is reasonable for 
low energy part of the spectrum, £_ < 112. In this case 
a slightly softer spectrum is produced at high energies 
for sufficiently high downstream temperature. This was 
shown in |l3| where such values of k and k\ have been 
employed. 



The solution F(v) in Eq.(Al) is in fact a one-sided 
(f|| < in the shock frame) isotropic (in the down- 
stream frame) part of the distribution function calculated 
at z = 0. To perform the matching with the high en- 
ergy standard (fully isotropic) power-law spectrum, a far 
downstream spectrum must be obtained since only the 
latter is isotropic in the downstream frame also at lower 
energies. The necessary formulae are given in ]l3| |. For 
the purpose of comparison with the simulation spectra 
given in p3| , this far downstream spectrum is pitch-angle 
averaged in the shock frame and drawn in Fig. M. 
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FIG. 2. Phase space of particles moving in the wave behind 
a shock. The shaded region corresponds to the adiabatically 
escaping particles. 
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Pitch angle scattering 
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FIG. 1. Velocity space of particles just in front of the 
shock in the shock frame of reference. The "beam" parti- 
cles schematically represent the leakage from the thermal dis- 
tribution downstream. After being scattered in pitch angle 
first upstream and then downstream, as shown by the ar- 
rows, these particles appear in front of the shock again but 
at higher energies and overlapping with beam particles. They 
are shown as "beam 2". 




FIG. 3. Contourplot of A for v = 0.9 and e = 0.2. 
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FIG. 6. The actual injection rate calculated from Eq. J7c| ) 
(solid curve); the approximations given by Eq. (|73| ) (dotted 
curve) and Eq. ^7% (dashed curve). 




FIG. 5. Injection efficiencies of different species normalized 
to proton efficiency as functions of mass to charge ratio and 
for different wave amplitudes. 
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FIG. 8. Particle spectra behind the shock, pitch angle av- 
eraged in the shock frame. The squares are from the hybrid 
simulations p3| . Thin line is a correspondent Maxwellian fit 
which is taken as a source term /m in Eq.(^). The solution 
of this equation is shown with the heavy line. The dashed 
line shows the solution of the same equation for t = 1. The 
dotted-dashed line indicates the slope of the spectrum appro- 
priate for high energy particles and shock compression of 4. 
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